Theory for Gossamer and Resonating Valence Bond Superconductivity 



m ■ 

o : 
o . 

(N ■ 

<: 
o : 

(N ■ 



i 

CO 

-4-* 
c3 



i 

c 

o 
o 



> 

00 
OS 

m 
oo 

o 
m 
O 

c3 



S3 
O 

o 



X 
S3 



J. Y. Gail 1 , F. C. Zhang 2 ' 3 , Z. B. Su 1 

1 Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, China 

2 Department of Physics, The University of Hong Kong, Hong Kong 
3 Department of Physics, University of Cincinnati, Cincinnati, Ohio 45221 Q/ 

We use an effective Hamiltonian for two-dimensional Hubbard model including an antiferromag- 
netic spin-spin coupling term to study recently proposed gossamer superconductivity. We formulate 
a renormalized mean field theory to approximately take into account the strong correlation effect in 
the partially projected Gutzwiller wavefucntions. At the half filled, there is a first order phase tran- 
sition to separate a Mott insulator at large Coulomb repulsion U from a gossamer superconductor 
at small U. Away from the half filled, the Mott insulator is evolved into an resonating valence bond 
state, which is adiabatically connected to the gossamer superconductor. 



1. Introduction 

Since the discovery of high temperature superconduc- 
tivity in the cuprates 0, there have been a lot of 
theoretical efforts trying to understand the microscopic 
mechanism for it. One of the scenarios was initiated by 
Anderson who proposed the idea of resonating va- 
lence bond (RVB) state for the observed unusual proper- 
ties in these compounds. A minimum model for cuprates 
was argued to be 2-dimensional Hubbard or its equiva- 
lent t ~ J model in the large U limit HQ. In the RVB 
picture, each lattice site is either unoccupied or singly oc- 
cupied by a spin-up or spin-down electron. The spins are 
coupled antiferromagnetically without long range order. 
The charge carriers move in the spin background and con- 
dense to a superconducting state 0, 13 El 13 EJ E3 • 
In this scenario, the undoped cuprate with density of 
one electron per site is a Mott insulator, and the super- 
conductor is viewed as a doped Mott insulator. Many 
experimentally observed properties in cuprates, such as 
the d-wave symmetry in superconductivity [lii . Il ll| ,the 
pseudogap phenomena 15], and the linear doping de- 
pendence of the superfluid density in the underdoped re- 
gion E3 , seem to be consistent with the RVB mean field 
theory. On the other hand, while mean field theories 
and variational calculations show the superconductivity 
in the doped Hubbard or t — J models, more direct nu- 
merical calculations on these models remain controversial 
and have been unable t o p rovide unambiguous answers 
to this question [HI EI El El El E3 

Very recently, Laughlin has proposed an interesting 
new notion, the gossamer superconductivity, for high T c 
superconducting Cu-oxides |23j|. In a gossamer super- 
conductor, the superfluid density is tenuous, in contrast 
to the conventional superconductor. He proposed an ex- 
plicit many body wavefunction for that state, which is a 
partially Gutzwiller projected BCS state (Eq. (2) below). 
The partial projection operator enables one to construct 
its inverse operator. Using these operators, Laughlin has 
further proposed a Hamiltonian, for which the partially 
projected BCS state is an exact ground state. Exact solu- 



tions play important role in many physical problems. By 
expanding that Hamiltonian, Laughlin showed that the 
superconducting ground state requires a large attractive 
interaction in addition to a large on-site Coulomb repul- 
sion. This raises the question if the Hubbard or t — J 
models capture the basic physics in cuprates. 

In a previous paper, one of us argued that the effec- 
tive Hamiltonian of the Hubbard model acting on the 
Gutzwiller's wavefunction should include a spin-spin cou- 
pling term, and the on-site Coulomb repulsion plays both 
the roles in projecting out (partially or completely) the 
double occupied state and in generating an attractive 
pairing interaction |24j. Although any variational calcu- 
lation cannot make a conclusion about the exact ground 
state, it is clear that the spin-spin coupling generated 
from the Hubbard U should capture some basic features 
in cuprates including its superconductivity. Nevertheless, 
Laughlin's idea of gossamer superconductivity is interest- 
ing. The wavefunction he proposed enables us to study 
the phase transition between a Mott insulator and a gos- 
samer superconductor at the half filled electron density 
and to study the strong coupling RVB state from a new 
viewpoint: namely its adiabatic continuation to the in- 
termediate coupling gossamer state. 

In this paper, we use an effective Hamiltonian (Eq. (1) 
below) for 2-dimensional Hubbard in square lattice to 
systematically study the partially projected Gutzwiller 
wavefunction. We are interested in the competition be- 
tween the Mott insulator and the superconductor. Here 
we shall neglect the possible antiferromagnetism in the 
model, which will be a subject in a future publication. 
We use Gutzwiller's approximation to replace the strong 
correlation in the projection by a set of renormalized fac- 
tors, and to use a renormalized mean field theory to study 
the ground state and the elementary excited states of 
the system. Our main results can be summarized below. 
At the half filling, the ground state is a Mott insula- 
tor at large U, and a gossamer superconductor at small 
U. The transition is first type in the physically inter- 
esting parameter region. The charge carrier density and 
the superconducting order parameter change discontinu- 
ously from zero in the Mott insulating phase to a finite 
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FIG. 1: Snap shots for a) A spin liquid Mott insulator; b). 
An RVB superconducting state; and c). A gossamer supercon- 
ducting state. Each bond connects an up-spin and a down-spin 
and open circle respresents the hole. 

value at the critical value of U. Away from the half filled, 
the gossamer superconducting state changes continuously 
from its state at the half filled, while the Mott insulat- 
ing phase becomes RVB superconducting. The gossamer 
and RVB suerconducting states have the same pairing 
symmetry, and their superconducting order parameters 
are both suppressed by a unified renormalized factor, 
which quantitatively characterizes the smallness of the 
superfluid density. Therefore, the gossamer and RVB su- 
perconducting states are adiabatically connected. The 
gossamer superconducting state at the half filled may be 
viewed as a RVB state with equal number of indepen- 
dent empty and doubly occupied sites. These empty and 
doubly occupied sites provide "parking space" for singly 
occupied electrons to move through the lattice. From 
this point of view, the relative reduction of U respected 
to the kinetic energy, which may be realized by applying 
the pressure (2j| , plays a similar role as the chemical dop- 
ing. The gossamer superconductivity may have already 
been realized in organic superconductors 26]. In Fig. 1, 
we schematically show snap shots for a Mott insulator, 
an RVB state, and a gossamer superconducting state. 

This paper is organized as follows. In Section 2, we in- 
troduce the model and the variational wavefunction. In 
Section 3, we formulate a renormalized mean field the- 
ory to study the variational wavefunction. Section 4 is 
devoted to the phase transition between the Mott insu- 
lator and the gossamer superconductor at the half filled. 
Detailed discussions on the gossamer and RVB supercon- 
ductivity are given in Section 5. The paper is concluded 
with a summary in Section 6. 



2. The Model and the Variational Wavefunction 

We study an effective Hubbard Hamiltonian in a square 
lattice, 

H = Ht+Hs+Hu (I) 
H t = - (*« c L c jV + h.c.) 

<ij><J 

H s = J Si ■ Sj 

<ij> 



i 

In the above equations, Ci a is the annihilation operator of 
an electron of spin a at the lattice site i, and n ia = cj^Cjo-. 
The sum is over the nearest neighbor pairs of < ij >, and 
U > is the intra-site Coulomb repulsion. Without loss 
of generality, we consider the case t > 0. In this Hamilto- 
nian, we have introduced an antifcrromagnctic spin-spin 
coupling term to account for the effect of the virtual elec- 
tron hopping process. In the large U limit, J as At 2 /U . 
This model may be viewed as an effective Hamiltonian 
of the Hubbard model. The inclusion of the antiferro- 
magnetic spin coupling appears consistent with the weak 
coupling renormalization group analysis |27| . and is ap- 
propriate in the variational approach studied here. In the 
limit U — > oo, the model is reduced to the t — J model. 
Very recently, a similar form of the Hamiltonian has 
been derived by using two subsequent canonical trans- 
formations starting from the Hubbard model at large U 
limit |2^. In the Hilbert space with the fixed number of 
the doubly occupied electron sites on the lattice, it has 
been shown that the Hubbard model may be mapped 
onto Eq. (1) with the constraint that the hopping pro- 
cess in Ht is limited to the corresponding Hilbert space, 
namely it does not change the electron double occupa- 
tion. Here we shall consider Eq.(f) from a phenomeno- 
logical point of view, and study its solutions within the 
framework of Gutzwillcr's variational approach. We shall 
consider J to be an independent parameter. 

Due to the perfect nesting and the van Hove singularity 
in the density of state, the ground state of Hamiltonian 
(I) at the half filling (electron density n = 1 per site) is an 
antifcrromagnet for arbitrarily small value of U even in 
the absence of the spin-spin coupling term. The spin-spin 
coupling further enhances the magnetism. In this paper, 
however, we shall focus on the insulating and metallic 
phases of the problem, and will not include the magnetic 
long range order. 

We study the model using a variational trial wavefunc- 
tion proposed by Laughlin 23], 

|*GS> = n^scs) (2) 

n« = - an^riii) ( 3 ) 

i 

with \^>bcs) a BCS-type superconducting state, given by 

i* S c5>-nK^+^4 T c U )|o> - (4) 
k 

where |0) is the vacuum, and and are variational 
parameters, satisfying the condition 

w\ 2 + kl 2 = i- 

n Q is a projection operator to partially project out the 
doubly occupied electron states on each lattice site i. 
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The state |\&gs) may be considered as a generalization of 
the previously studied partially projected non-interacting 
electron state [2|| yfj, |3lj to include the superconducting 
state. In the limiting case «£t>£ = 0, \^bcs) is reduced 
to the non-interacting electron state, and 

I*gs> -> n Q |* Fi ), 

where I^^l) is the ground state of the non- interacting 
electron system, given by \&fl) = II ct C £ CT |0), with 

the product running over all the fc's within the Fermi 
surface. |\&gs) is a natural generalization of the usual 
BCS state to strongly correlated systems. It connects 
the usual BCS state to the RVB state, characterized by 
the parameter a, which takes the value between and 1. 
a = is a normal BCS state. At a = 1, the projection 
operator projects out all the doubly occupied electron 
states, and I^gs) is reduced to the RVB state At 
the half filling and at a = 1, each lattice site is occupied 
by a single electron, and the system is a Mott insulator. 
Therefore, the wavefunction |\&gs) can be used to study 
superconductor-insulator transition. 



3. The Renormalized Mean Field Theory 

We now proceed the variational calculations to deter- 
mine the parameters a and uz, wg. Without loss of gen- 
erality, we consider the electron density n < 1. The 
variational energy per site E = (H) is given by 



E = Ud+ (H t ) 



(H.j 



(5) 



where d = (n^riii) is the average electron double oc- 
cupation number. (Q) is the expectation value of the 
operator Q in the state |^gs)- F° r briefness, (H t ), (Hj) 
stand for their average values per site, d is a function of 
a, and < d < 1/4. The first term in Eq. (5) is the 
intra-site Coulomb interaction energy, while the second 
and the third terms are the average kinetic and spin-spin 
correlation energies, respectively. 

The variational calculations can be carried out using 
variational Monte Carlo method H HI HI Ei- Here 
we use the renormalized Hamiltonian approach to treat 
the projection operator approximately [7j. In this ap- 
proach, the effect of the projection operator is taken 
into account by a classical statistical weighting factor, 
which multiplies the quantum coherent result of the non- 
projected state. This method (Gutzwiller method here- 
after) was first proposed by Gutzwiller [3(J , and has been 
applied to study strongly correlated systems by many 
others 

HE El- Let (Q) o be the expectation value 
of Q in the state \^bcs), then the hopping energy and 
the spin-spin corerlation in the state I^gs) are related 
to those in the state \^bcs) by, 




FIG. 2: The Gutzwiller's renormalization factors, g t and g s , 
as functions of the double occupation number d, obtained 
from Eqns. (6). S = 1 — n e . 



(Si ■ Sj) — g s (Si ■ Sj)o 



(6) 



The renormalized factors gt and g s are determined by 
the ratio of the probability of the physical processes in 
the states |\&gs) and \^bcs)- Following the counting 
method described in the literature [?J, we have 



gi 

9s 



(n - 2d)(Vd+ Vl-n + d) 2 
(1 -n/2)n 
(n - 2d) 2 



(l-n/2) 2 n 2 



(7) 



The expression for g t is the same as in the early liter- 
ature 29]. In the limit d — 0, Eqns. (7) recover the 
previous results derived for the t — J model [jj • These 
renormalized factors quantitatively describe the correla- 
tion effect of the on-site repulsion. g t < 1, and g t << 1 
at small d and small (5, representing the reduction of the 
kinetic energy due to the projection. 4 > g s > 1, and 
g s = 4 at d = and 5 = 0, representing the enhance- 
ment of the spin-spin correlation due to the projection. 
In Fig. 2, we plot gt and g s as functions of the double 
occupation number d for various electron densities. 

In terms of these renormalization factors, we can define 
a renormalized Hamiltonian given by 



H' = g t H t + g s H s + H L 



(8) 



The expectation value of H in the state I^gs) can be 
evaluated in terms of the expectation value of H' in the 
state Y$bcs)- We obtain, 



E = (H') = Ud + g t (H t )o + g s (Hj) 



(9) 



(cLc,- 



9t(4 a c ja ) 



In the renormalized Hamiltonian approach, the origi- 
nal variational parameters {a,?^,?^} are transformed 
into the variational parameters {d,v^u^\. There is one 
to one correspondence between a and d. Within the 
Gutzwiller approximation, one can analytically calculate 
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d = (Tiiiriii), and one finds |3Cl l3l|. 
(1-a) 



2 d(l — n + d) 



(n/2-d) 2 



(10) 



The energy in Eq. (9) is evaluated for the fixed number 
of electrons N e . We introduce a Lagrangian multiplier /i, 
and define 



K = H' - fl(J2n ia - N e ). 



(11) 



We then have E = (K)o, being subject to the condition 
d(K) /dfl = 0, or 



(12) 



Below we consider the case and to be real. Evalu- 
ating Eq. (11), we obtain (lattice constant =1), 

E = Ud+fL + 2j2(9te%-ij>)v 2 n 

k 

+ Yl V k-9 ( v k 2v £> 2 + u k v k u k> V P ) (I 3 ) 
%,P 

where 



V k = - -^9 s J (cos k x + cos k y ) 
Ej: = — 2i(cos k x + cos k y ) 



(14) 



Carrying out the variational procedure with respect to 
M£ and ug, we obtain 



where 



(15) 



(16) 



The variational parameters and xg are related to the 
particle-particle amplitude and the particle-hole par- 
ing amplitudes xg by, 

Aj: = A x cos k x + A y cos k y 

Xfe = eg- {xx cos k x + x y cos k y ). (17) 

In the above equations, we have introduced two correla- 
tion functions in the unprojected state \^BCs)y 

At = (CiJ.Cj+ T ,! — Ci|C 4+Tj |)o 
a" 

with r = x, y, the unit vectors on the lattice, and 

eg = [-2g t t(cosk x + cosk y ) - fi]/(3g s J/A). (19) 



For the <i-wave pairing state, which has the lowest energy 
within this class of states as suggested in the previous 
studies for the t — J model 0, Ifl lift Infil l , we have A x — 
—A y = A, and Xx — Xy — X- A and x are determined 
by the coupled gap equations, 



A = ^(cosfe x )Ag/£g 
k 

X = -J2( cosk ^)Xk/ E k 



(20) 



These gap equations must be solved simultaneously with 
the hole concentration equation, Eq. (12), which can be 
rewritten as S — ^ xg/^gi with S = l—n e . The variation 

fc 

with respect to d leads to the equation 

In terms of x and A, the energy is given by 

E = Ud-4ng t t X -(3g s J/4)(A 2 + x 2 ) (22) 

where x and A are the solutions of the gap equations, and 
both are functions of d. In the case there are multiple 
solutions for d from Eq. (21), the ground state is deter- 
mined by the global energy minimum. Alternatively, we 
may solve the gap equations for given values of d, and 
calculate E(d) to find the optimal value of d to deter- 
mine the ground state and the ground state energy. The 
chemical potential of the system, \l = dE/dn, is given by 



dgt 
dn 



to 



dg s 
dn 



H. 



JO 



(23) 



Note that chemical potential here is different from the 
Lagrangian multiplier jl in the renormalized mean field 
theory. This is because the renormalized factors g t , g s to 
be also functions of electron density n. 



4. Mott Insulator-Gossamer Superconductor 
Transition 

In this section, we discuss the variational solutions at 
the half filled case. At the half filling, the trial wave- 
function |\&gs) describes either a Mott insulator if a = 1 
(i.e. d = 0), or a superconducting state if a < 1 (i.e. 
d > 0). If a is close to 1, or d is very close to zero, |^gs) 
describes a gossamer superconducting state. 

We expect a Mott insulator at large U and a supercon- 
ducting state at small U. This can be examined qualita- 
tively without carrying out the quantitative calculations. 
At the half filling, g t = 8(1 - 2d)d, and g s = 4(1 - 2d) 2 . 
Eq. (21) becomes 

17 + 8(1 - 4d){H t ) - 16(1 - 2d)(Hj) = 0. (24) 
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Since both (H t )o and (-Hj)o are finite, there will be no 
solution of Eq. (24) if U is sufficiently large. This indi- 
cates that the ground state corresponds either d = or 
d = d max , the allowed maximum value of d. The repul- 
sive nature of U excludes the latter, and it follows that 
the Mott insulating state with d = is the ground state. 
We believe that the qualitative result for the existence of 
the Mott insulating phase at large but finite U is robust. 
Note that in the Gutzwiller's wavefunction, the doubly 
occupied site and the empty site are not correlated. At 
the half filling, d represents the carrier density n* and 
is proportional to the Drude weight in the a.c. conduc- 
tivity, n*e 2 /m*, with m* the effective mass. We remark 
that the parameter d in our Gutzwiller approach is dif- 
ferent from the usual double occupation number d (for 
example, the double occupation calculated in the exact 
diagonalization of a finite size system). In the latter case, 
d also includes the contribution from the virtual hopping 
process. Therefore, the double occupied site is bound to 
the empty site, and the double occupation number d does 
not represent the mobile carriers. 

In the insulating phase, the ground state is the same 
as that of the Heisenberg model. Within our theory, the 
ground state energy is given by 

E Q = -3J(A 2 ) + X 2 ) (25) 

with 

A = xo = C/V2 

= i ^cos 2 k x + cos 2 k y = 0.339. (26) 

At small U, one expects a metallic ground state, except 
in the special cases due to the band effect, such as having 
the Von Hove singularity and perfect nesting in H t . In 
this paper, we will not consider the special band effect. 
Therefore, we expect a metal-insulator transition at a 
finite U = U c in the general case, with the metallic phase 
to be superconducting provided that ugug ^ 0. This is 
the Mott insulator-gossamer superconductor transition. 

We now discuss the phase transition in details. We 
solve the gap equations for the fixed d and determine 
the transition point U c and the nature of the transition. 
The phase diagram in the parameter space U and J/t is 
plotted in Fig. 3. The critical U c separates the Mott in- 
sulating phase from the gossamer superconducting phase. 
We can choose the mobile carrier density as the order pa- 
rameter, which is proportional to d. The phase transition 
is classified as the second type if d — > and the first type 
if d — > d c > as U — ► U c within the metallic phase. This 
classification is consistent with the usual zero tempera- 
ture quantum phase transition, where the nature of the 
phase transition depends on the continuity or disconti- 
nuity of the order parameter. We find the transition to 
be the first order at < J/t < r) c , and the second order 



Uc/t 




12 3 



J/t 

FIG. 3: Phase diagram at the half filling. 

at r/c < J/t, with r/ c w 2. At J — 0, the present theory is 
reduced to the Brinkman-Rice theory for metal-insulator 
transition | 29| for the projected non-interacting electron 
state. In that case, we find U c /t — 128/7T 2 . From Fig. 3, 
we see that U C (J — > 0) = U C (J = 0), so that the critical 
value of U is continuous at J — 0. However, the transi- 
tion is the second type at J = 0, while it is the first type 
for any small but finite J/t. 

Let us first discuss the first order phase transition in 
the region < J/t < 2. In Fig. 4, we show the energy 
E as a function of d for several values of U at a typical 
parameter J/t — 1/3. E is not a monotonic function 
of d for U near U c . As d increases from d = 0, E first 
increases linearly, then decreases, then increases again. 
There is a local energy minimum around d = 0.02, which 
develops and becomes a global minimum as U approaches 
U c from the insulator side. The local minimum E(d c ) at 
d = d c represents a metallic solution, and E(d = 0) rep- 
resents an insulator solution. The critical value for the 
Mott insulator and gossamer superconductor transition 
is determined by the condition E{d c ) — E{d — 0). From 
Fig. 4, we have U c /t = 10.23 for J/t =1/3. At U > U c , 
d = 0, and the ground state is an insulator. At U < U c , 
d> d c fa 0.02, and the ground state is a gossamer super- 
conducting state. In Fig. 5, we plot d as a function of 
U. d is approximately linear in U till the transition point 
U c . The discontinuity in d is about 0.02. We conclude 
that the Mott insulator-gossamer superconductor phase 
transition in this most relevant region is first type. The 
carrier density is discontinuous at the phase transition 
point. Since d is proportional to the carrier density, this 
type of first order transition should be observable in the 
electric transport or in the a.c. conductivity measure- 
ments. 

For large ratio of J/t, our calculations show that the 
phase transition is second order. This is illustrated in 
Fig. 6 for E v.s. d in the case of J/t = 3. The transition 
occurs at U c = 2.58i, and d changes continuously across 

Uc- 

A special case is J = 0. In this limit, |^gs) = 
IIcl^Fi}) a nd our theory is reduced to the previous one 
for the projected Fermi liquid state. The energy in Eq. 
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FIG. 4: The energy E as a function of d for U around the 
critical value U c at S = 0. The ratio J/t = 1/3. 



0.25 




5 10 15 

U/t 

FIG. 5: The double occupation number d as a function of U 
at 5 = and J/t=l/3. 

(9) at the half filling becomes 

i? = [/rf — 2g t | cos fc-r + cos k y \ 

k 

= Ud-128d{l-2d)t/n 2 (27) 

From this we find U c /t = 128/7T 2 w 13. d is continuous 
at U c so that the transition is second type. 

Our result on the first order phase transition in the 
physically interesting region (small but non-zero J/t) is 
somewhat unexpected. We argue that the first order 



E(d)/t 



-2.045 
-2.05 


U/t=2.4(5=0.0) 
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FIG. 6: The energy E as a function of the double occupation 
number d for several values of U around the critical value U c 
with 8=0 and J/t=3. 



transition between the Mott insulator and the gossamer 
superconductor is due to the interplay of the kinetic and 
spin-spin correlation energies. This interplay was not in- 
cluded in the previous study of the Gutzwiller approach 
but is taken into account here. To illustrate the effect 
of the interplay to the nature of the phase transition, we 
consider the limiting case < J/t << 1, and expand the 
energy E of Eq. (22) at n = 1 for small d, 

E{d) =E + {U - U c0 )d - 13d 2 + 0(d 3 ) (28) 

where E is the energy at d = given by Eq. (25), U c0 — 
l6V2Ct - 12C 2 J is the solution of dE/dd = at d = 0, 
given by Eq. (24). j3 = [32(dx/dd) ld=g - 32V2C)]/j. The 
J-dependence in (3 has been neglected since J/t << 1. 
Note that the kinetic energy is proportional to x- As d 
increases from 0, x tends to increases from xo = C/y/2 to 
gain more kinetic energy. Therefore, dx/dd > 0. In the 
limit J/t « 1, weh&ve dx/dd = (x~Xo)/d cat/ J » 1, 
hence the first term in the expression for (3 dominates 
and j3 > 0. This demonstrates that d — is a local 
maximum in energy at U = U c q, and the phase transition 
occurs at a large value of U corresponding to d > as 
numerically shown in Fig. 4, hence it is a first order 
transition. Numerically, we find that (3 — 34.8472, for 
J/t = 1/3. 

It is interesting to compare the gossamer superconduc- 
tor - Mott insulator transition with the metal-insulator 
transition studied in previous literature j^. In the 
Brinkman-Rice theory, the transition is second order. 
In that theory, the system approaches the insulating 
phase, the effective mass m* — > oo. In the gossamer 
superconductor-Mott insulator transition with small ra- 
tio of J/t, the insulating phase is not characterized by 
the divergence of the effective mass. We estimate the ra- 
tio of the effective mass to the band mass (1/t) at the 
metallic side of transition point to be 1/gt ~ l/(8d) ~ 6. 

The first order phase transition between metal and in- 
sulator was first pointed out by Peierls [13] and by Lan- 
dau and Zeldovich j^, and examined in more great de- 
tail by Mott 1^3. In their theory, an electron is always 
bound to a positive charge due to the long range Coulomb 
attraction, and the transition of a metal to an insula- 
tor at zero or very low temperatures occurs at a finite 
critical electron density, and must be first type. It is 
interesting to note that the on-site repulsion also leads 
to the first order transition between a specific type of 
metal (superconductor) and an insulator studied in the 
present paper, where the long range Coulomb force is not 
included. We also note that Florencio and Chao 0] in- 
vestigated the metal-insulator transition of the Hubbard 
model using Gutzwiller's wavefunction by including anti- 
ferromagnetism and found the transition to be first type. 
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FIG. 7: The double occupation number d as a function of 5 
at the ratio J/t — 1/3 for several Vs. 



5. Gossamer and RVB superconductivity 

In this section, we discuss the superconducting state 
at the half filled as well as away from the half filled. 
Note that at 6 > 0, |^gs) always describes a metallic 
state. To make the terminology clearly, we shall call the 
superconducting state at U < U c to be the gossamer 
superconductor [2^, |24j, and the doped Mott insulator 
(U > U c and 8 > 0) to be the RVB state 0. 

We begin with the discussion of the double occupation 
number d as a function of the hole concentration. Wc 
solve the gap equations and find the optimal value of d. 
The results are plotted in Fig. 7. We find that d is al- 
ways non-zero at 5 > 0, even in the region U > U c . This 
suggests that the doped Mott insulator is described by 
a partially projected state (a < 1 in Eq. (2)). Nev- 
ertheless, d is very small for U 3> U c . As we can see 
from Fig. 7, d varies from to 0.01 for U/t = 15, which 
corresponds to U/U c w 1.5. The non-zero value of d at 
5 > may be understood from the variational equation 
(21), which determines d. At 6 > 0, dg t /dd\ d=0 — > oo. 
Therefore, d = cannot be a solution of the equation, 
and d must be finite. It remains to be seen if this result 
is due to the Gutzwillcr's approximation used in our cal- 
culation. It will be interesting to further examine this 
issue using other methods such as the variational Monte 
Carlo method. 

From Fig. 7, we also see that as 6 increases from 0, 
d increases for large U while d decreases for small U. 
The latter may be understood as follows. In the small 
U case, the correlation becomes less important, and the 
qualitative feature between d and 5 becomes similar to 
the uncorrelated state. For the uncorrelated Fermi liquid 
state, d = (1 — S) 2 /4, so that d monotonically decreases 
as 5 increases. 

While d is a smooth function of 6 for most values of 
U in our study, there is a narrow region in U above U c , 
d changes discontinuously at a very small S. In Fig. 8, 
we show the energy E v.s. d for U = 10.235t, slightly 
above U c = 10.23t, for four values of S. At 6 = 0, d = 



m 




-0.237 



. 0.01 0.02 0.03 0.01 0.02 0.03. 




FIG. 8: Energy E as a function of d for U = 10.235i > U c at 
the ratio J/t = 1/3, for several values of S. 



corresponds to the global energy minimum, while there 
is a local minimum around d — 0.02. As i5 gradually in- 
creases, the positions of the two minima change smoothly 
and their corresponding energies reverse their order. In 
this region, the optimal value of d jumps. This region is 
found very narrow: 10.23t < U < 10.235i, however. 

We now discuss the superconducting order parameter. 
The superconducting order parameter of the state I^gs) 
is defined by, for the d-wave pairing, 



A sc (t) = (c li c 4+rT ) - (c lT c i+ri ), 



(29) 



and A sc = A sc (x) = — A sc (y). We shall adopt the 
Gutzwiller approximation to calculate this quantity. In 
analogy to the derivation for the hopping energy in 
Eq.(6), we find that 



< Cij.c l+rT > = g t < c a c i+rT > 



(30) 



Therefore the order parameter A sc is related to the vari- 
ational parameter A in the gap equations by 



St A 



(31) 



In Fig. 9, we show our results for A sc and A as functions 
of S for three values of U: well above U c , at U c , and well 
below U c - Note that at U = 15t >> U c , A sc = at 
5 = 0, although A takes a maximum. This is consistent 
with the Mott insulating ground state. At 5 increases, the 
kinetic energy plays more important role in comparison 
with the spin-spin correlation energy, and A decreases 
monotopically. However, A sc shows a non-monotonic 
dome shape, and it first increases to reach a peak be- 
fore it drops for larger U. Also note that at U c the Mott 
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FIG. 9: Variational parameter A and superconducting order 
parameter A sc as functions of the hole concentration 5 for 
three values of U at ratio J/t=l/3. 



insulator and gossamer superconducting state are degen- 
erate at 8 — 0, and the gossamer superconducting phase 
continuously evolves into the metallic phase at 8 > 0. 
Shown in the figure for U = 10.23£ « U c is the metallic 
phase. The non-zero value of A sc at 8 = indicates the 
transition to be the first order. 

In Fig. 10, we show A sc and A as a function of U for 
<5 = 0. From this figure we can see that for U < U c , 
A sc and A both increase with U . At U > U c A takes a 
maximum while A sc = 0. Note that at small values of U, 
the effective Hamiltonian does not represent the original 
Hubbard model, and one should be cautious to interpret 
the results at small U. 

It is interesting to point out that the superconducting 
order parameter A sc in both the gossamer and RVB su- 
perconducting states are characterized by the variational 
parameter A and a small renormalized factor g t . They 
have the same pairing symmetry, and the two states are 
adiabatically connected to each other. From this point 
of view, the gossamer and RVB superconductors are the 
same. In Fig. 11, we present a schematical diagram for 
the Mott insulator, gossamer and RVB superconducting 
states. In the parameter space of Coulomb repulsion U 
and of the hole concentration 8, there is a line at 8 = 
and at U > U c for the Mott insulating phase. In the 
region U > U c and 8 > 0, it is the RVB superconducting 
phase as Anderson proposed {4j. In the region U < U c , it 
is the gossamer superconducting phase in which both the 
half filled and non-half filled states are superconducting. 

While the gossamer and RVB superconducting states 
are essentially the same, the chemical potential /j, in the 
gossamer superconducting state is continuous at 8 = 



FIG. 10: Variational parameter A and superconducting order 
parameter A sc as functions of U at S — and J/t = 1/3. 
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FIG. 11: Schematical phase diagram for the Hamiltonian in 
parameter space U and 5. 

because of the metallic phase, while /i is discontinuous 
at 8 = because the state at 8 — is an insulator and 
the state at any small but finite 8 is a metal within the 
present theory. 

Below we shall study fi quantitatively. At 8 = 0, /i = 
U/2 by electron-hole symmetry. At other value of 8, we 
calculate /i using Eq. (23) after solving the gap equations. 
In Fig. 12, we show /i as a function of 8. As we can see 
from the figure, /x = U/2 at 8 = 0, and is continuous 
for U < U c = 10.23i. There is a discontinuity in fi for 
U > U c at 8 = 0. At U > U C) the chemical potential is 
shifted from U/2 at the half filled to the lower Hubbard 
band away from the half filled. To see this more explicitly, 
we define A/x = n(d — > 0) — fi(d = 0). Afj, as a function 
of U is plotted in Fig. 13. As U decreases, A/x decreases 
monotonically and reaches a finite value at U = U c + + , 
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feature was first found for the large-U limit Hubbard 
model 0, and is consistent with the "high energy pseu- 
dogap" seen in the angular resolved photoemission exper- 
iments and the experimentally observed superconducting 
energy gap [4l|, |42> |43| . Here we show that this feature 
appears also in the gossamer superconductor. 



6. Summary 



FIG. 12: The chemical potential /i as a function of the hole 
concentration 8 for several values of intra-site Coulomb repul- 
sion U with J/t = 1/3. 



Au/i 




8 9 10 11 12 13 14 15 



u/t 

FIG. 13: The discontinuous in the chemical potential as a 
function of intra-site Coulomb repulsion U for J/t=l/3. 

then drops to zero at U = U c — + . The discontinuity of 
A/i at U c is related to the first order phase transition. 

Finally, we briefly discuss the excited states. In the 
context of the Gutzwiller trial wavcfunction, the excited 
states were discussed by Zhang et al. Q for the RVB 
state, and recently by Laughlin for a gossamer super- 
conductor Harniltonian|23| . Following Zhang et ah, the 
quasiparticle states can be defined by 

The quasiparticle energy Ep is defined to be the dif- 
ference of the expectation values of K in Eq. (32) in 
this state and in the ground state |^gs)- We use the 
Gutzwiller method to calculate the energy and obtain 

Ep=(3g s J/4)^ X }+^ (33) 

At the vectors p satisfying Xp- = (Eq. (17)), we have 
Ep = A | cos k x — cos k y \ . Therefore, the quasipartical en- 
ergy is proportional to the parameter A, and not renor- 
malized by the factor g t , which is very different from 
the superconducting order parameter. From Fig. 9, 
our theory predicts the quasipartical energy to be max- 
imum at S = 0, and decreases as doping increases. This 



We have used the Gutzwiller variational method to 
study an effective Hamiltonian for 2-dimension Hubbard 
model. Based on Gutzwiller approximation, we have dis- 
cussed the case both at the half filled and away from 
the half filled. At the half filled, there is a first order 
phase transition to separate a Mott insulator at large 
Coulomb repulsion U from a gossamer superconductor 
at small U.This is very interesting. It suggests that the 
on-site Coulomb repulsion can lead to the first transition 
between a specific type of metal and an insulator. The 
double occupation number d which is proportional to the 
carrier density changes discontinuously from zero in the 
Mott insulator phase to a finite value at the phase tran- 
sition point([/ = U c )- So we expect that this type of 
first order transition should be observable in the elec- 
tric transport or in the a.c. conductivity measurements. 
Away from the half filled, the Gutzwiller variational state 
is always metallic. The gossamer superconducting state 
changes continuously from its state at the half filled, 
while the Mott insulating phase becomes RVB supercon- 
ducting. The gossamer superconductor is similar to the 
RVB suerconducting states with the same paring sym- 
metry, and showing the pseudogap. Their major differ- 
ence is on the position of their chemical potential. The 
Gutzwiller method we used in this paper has previously 
been tested in good agreement with variational Monte 
Carlo method pj 0- We believe that the qualitative 
conclusions obtained here should be reliable, and refined 
numerical calculations such as variational Monte Carlo 
calculations will be interesting to examine the problem. 

There are other questions that require further investi- 
gation such as the effect of the antiferromagnetism which 
is most plausible in this model. 

We wish to thank R. Laughlin, T. M. Rice, P. W. An- 
derson, Y. Yamashita, and Y. Yu for many interesting 
discussions. This work was partially supported by the 
US NSF grant 0113574, and by Chinese Academy of Sci- 
ences. 
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